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摘 E: 两 相 流 系统 中 ， 流 型 会 直接 影响 系统 的 摩擦 阻力 和 传 热 等 特性 ， 而 流 型 的 转变 过 程 会 受到 通道 的 形状 和 尺 
二 的 影响 。 常 温 常 压 下 ， 以 空气 和 水 为 介质 ， 对 八 种 不 同 尺寸 矩形 通道 中 竖 直 向 上 泡 状 流向 弹 状 流转 变 过 程 进行 了 
研究 。 结 果 表 明 ， 和 拖 形 通道 间隙 大 于 2 mm 时 ， 随 着 通道 水 力 直径 的 增 大 ， 泡 状 流向 弹 状 流转 换 边 界 左 移 。 在 窗 通 
道 范 围 内 ， 可 采用 Ishii (1977) 提出 的 基于 漂移 流 模 型 的 判定 准则 计算 流 型 转换 边界 ; 在 常规 矩形 通道 中 ， 转 换 准 
则 中 的 临界 空 泡 份 额 与 气泡 初始 尺寸 有 关 ， 选 用 赵 建 福 提 出 的 方法 计算 临界 空 泡 份额 ， 转 换 边 界 预 测 结果 与 实验 结 
果 符 合 程度 较 好 ; 在 利用 漂移 流 方法 计算 截面 含 气 率 时 ， 对 于 窦 矩形 通道 ， 分 布 系数 可 以 用 Ishii (1977) 提出 的 方 
法 计算 ， 而 对 于 常规 矩形 通道 ， 分 布 系数 为 1.2。 
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Criterion for Predicting Transition Boundary from Bubbly Flow to Slug 


Flow in Rectangular Channels 


JIN Guang-Yuan HAN Yue-Yang 
(School of Energy and Power Engineering, Northeast Electric Power University, Jilin 132012, China) 


Abstract: Flow pattern transition always affect the frictional resistance and heat transfer in two-phase 
systems, whereas it itself is also affected by the shape or size of the channels. The transition from bubbly 
flow to slug flow in eight different rectangular channels was investigated under normal condition. The 
results shows that the transition boundary from bubbly flow to slug flow moves to left as the diameter 
increasing when the channel gap is smaller than 2 mm. The transition criteria can be predicted from the 
model from Ishii (1977) in narrow rectangular channels. While for common channels, the critical void 
fraction in the transition criteria is relevant to the size of bubbles in entrance, for which calculation 
method of Jianfu Zhao (2005) is proposed, leading to the predicted transition criteria with a good 
agreement against the experimental result. When using of drift-flux model, the distribution parameter Co 
in channels is calculated with different forms, and the result shows that correction from Ishii (1977) can 
be used in narrow size of rectangular channels, while Co remains to be 1.2 in common channels. 


Key words: Flow transition; Bubbly flow; Slug flow; Void fraction; Initial size 


0 引言 直上 升 矩 形 通 道内 流 型 进行 研究 ， 发 现 通 道中 存在 
B 泡 状 流 、 弹 状 流 和 环 状 流 ， 当 直径 超过 10 mm 时 ， 

两 相 流 系统 中 ， 流 型 的 转换 会 直接 影响 流动 阻 管道 的 几何 形状 对 流 型 的 影响 很 小 ; Mishima 5?! 

力 和 换 热 特性 的 机 理 ， 基 于 流 型 的 经 验 公 式 或 半 经 通过 和 光照 相 方式 对 罕 矩 形 通 道中 的 流 型 进行 了 研 
验 公式 也 将 随 之 发 生变 化 ， 因 此 研究 各 种 通道 中 的 究 ， 在 间隙 为 1.0 mm 的 和 宗 窍 形 通道 中 搅 浑 流 没 有 
两 相 流 热 工 水 力 特 性 都 要 以 精确 地 流 型 判别 为 基 出 现 。Wilmarth 和 Ishii 等 人 BI 对 间隙 为 1 mm 和 2mm 
础 。 对 于 和 矩形 通道 中 流 型 的 研究 主要 在 国外 ， 见 文 的 水 平 罕 矩形 通道 ， 并 与 Ishii5I 和 Taitell 2x25 T- 18] 
献 [1-4]，Sadatomi 等 山 对 水 力 直径 超过 10 mm 的 竖 状 流 向 弹 状 流 过 渡 的 模型 进行 验证 ， 并 提出 新 的 用 
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于 计算 的 分 布 系数 Co. Hibiki 等 由 通过 对 间隙 为 
0.3—17 mm 的 矩形 通道 流 型 进行 研究 , 得 到 了 新 的 
关于 流 型 转换 的 准则 。 

对 于 和 矩形 通道 内 流 型 转变 准则 的 研究 已 开展 多 
年 ， 但 在 该 领域 得 到 的 结论 ， 尤 其 是 泡 状 流向 弹 状 
流转 变 的 准则 有 些 文献 给 出 的 结论 是 不 一 致 的 ， 本 
实验 选取 不 同 几 何 尺 寸 的 矩形 通道 ， 对 于 竖 直 向 上 
泡 状 流向 弹 状 流转 换 边 界 进行 研究 。 


] 实验 装置 


实验 在 常温 、 常 压 下 进行 ， 图 1 中 的 实 线 和 虚 
线 分 别 表示 水 回路 和 气 回路 。 实 验 段 分 别 为 截面 为 
43xl.6mm? (通道 一 ) . 43x3mm? (通道 二 ) 和 
43xl0mm? OGH3É —) 三 种 矩形 通道 ,， 水力 直 径 分 别 
为 3mm、5.6mm 和 16.2mm， 长 度 均 为 2000 mm, 
基于 Kandlikar 中 的 定义 ， 通 道 一 为 窄 矩 形 通 道 ， 通 
道 二 、 三 为 常规 矩形 通道 。 去 离子 水 由 离心 泵 提供 
驱动 力 ， 经 过 滤器 、 调 节 阀 和 混合 腔 进 入 实验 段 返 
回 水 箱 形 成 循环 ,空气 经 气泵 压缩 后 存在 储 气 负 中 ， 
使 用 时 经 过 减 压 阀 和 止 逆 阀 进入 混合 腔 、 实 验 段 后 
由 实验 段 上 方 汽水 分 离 装 置 进 入 大 气 。 混 合 腔 中 均 
布 若干 直径 为 1 mm 的 毛细 管 ， 确 保 实验 段 入 口 处 
气泡 大 小 接近 ， 发 泡 均匀 。 水 流量 和 气流 量 分 别 采 
用 不 同 质量 流量 计 测 得 ， 精 度 分 别 为 0.1 级 和 0.2 
级 。 距 实验 段 入 口 500 mm 及 1500 mm 处 安装 有 精 
度 为 0.2 级 的 压力 传感器 。 正 对 实验 段 两 测 压 点 中 
心 位 置 固定 安装 高 速 摄影 仪 ， 用 以 拍摄 气泡 行为 。 
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图 1 实验 系统 图 


Fig. 1 Schematic diagram of the experimental system 
2 泡 状 流 、 弹 状 流 流 型 的 判别 

本 实验 采用 高 速 摄影 方式 对 气 水 形态 进行 观 
测 ， 弥 补 人 眼 对 于 高 速 流体 辨识 不 清 的 缺点 。 文 献 
[1-4] 中 对 于 常规 矩形 通道 和 和 罕 矩 形 通 道中 的 泡 状 流 
和 弹 状 流 的 描述 是 不 同 的 ， 图 2 中 给 出 了 实验 中 不 
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同和 矩形 通道 中 泡 状 流 和 弹 状 流 的 示意 图 。 

在 常规 通道 中 ， 泡 状 流 内 气泡 是 以 彼此 分 离 的 
小 泡 形式 分 散在 液体 之 中 ， 在 较 小 间隙 的 通道 中 ， 
气泡 多 为 圆柱 形 ， 较 大 间隙 的 通道 则 为 球状 ， 弹 状 
流 的 气 弹 多 数 有 圆润 的 涉 部 和 扁平 的 主体 ， 气 弹 和 
液 相交 叉 分 布 ， 液 相 中 常 含有 大 量 弥散 的 小 气泡 。 
在 窄 矩 形 通 道中 ， 由 于 通道 间 际 的 限制 ， 泡 状 流 即 
使 在 气相 速度 较 低 时 也 难以 发 现 ， 气 泡 提 前 聚合 形 
成 帽 状 流 的 流 型 ， 气 泡 头 部 呈 半 圆 型 ， 下 部 近似 为 
水 平 状 ， 弹 状 流 与 常规 矩形 通道 近似 相同 ， 但 罕 通 
道 弹 状 流 的 液 相 部 分 内 一 般 仅 含 有 少量 的 气泡 ， 有 
时 气 弹 中 会 有 液 滴 出 现 。 


泡 状 流 弹 状 流 “ 泡 状 流 ER 
通道 通道 二 通道 三 
CD=3 mm) (D=5.6 mm) CD=16.2 mm) 


图 2 不 同和 矩形 通道 两 种 流 型 示意 图 


Fig.2 Two flow patterns in different kinds of rectangular 


channels 

3 泡 状 流 一 弹 状 流转 换 边 界 判定 方法 
3. 1 泡 状 流 一 弹 状 流转 换 边 弄 对 比 

选取 从 文献 2 和 文献 11 得 到 的 矩形 通道 中 泡 状 
流向 弹 状 流转 变 边 界 ， 并 与 实验 中 得 到 的 结果 进行 
比较 ， 如 图 3 所 示 (图 中 为 通道 间隙 ，w 为 宽 边 
长 度 , DD 为 当量 直径 ， 尺 寸 单位 均 为 mm)。 对 比 发 
现 和 矩形 通道 中 泡 状 流向 弹 状 流转 变 主要 与 通道 间 阶 
有 关 ， 当 通道 间隙 小 于 2 mm 时 ， 间 隙 为 0.3、1.0 
和 1.6 mm 的 转换 边界 几乎 重合 ， 当 通道 间隙 大 于 
2mm 时 , 随 着 通道 水 力 直径 的 增加 , 转换 边界 左 移 ， 
由 于 间隙 变 小 造成 的 束缚 减弱 ， 壁 面 的 剪 切 作用 也 
相应 变 小 ， 气 泡 间 更 容易 实现 碰撞 和 聚合 ， 弹 状 流 
出 现时 间 相 应 提前 。 
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Fig.3 Comparison of different corrections for distribution 
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parameter in rectangular channels 
3. 2 76:1 E SCA R F 3802 0] x TEE DU AY 
3. 2. 1 78 IA t SPI ce 158322 7| E E 
许多 学 者 都 采用 漂移 流 的 公式 作为 泡 状 流向 弹 
状 流 转变 的 判别 式 ，Jones 和 ZuberlS1 (1975) 提出 


基于 算 形 通道 的 转换 边界 判定 准则 ， 形 式 为 

1 1 
j=( -1) j, ——(0.23+0.13s / w) Apgw/ 
k= 5030 D ^c s/ Ww) Apgw! p, 


(D 

Ishii! (1977) 得 到 了 和 矩形 通道 泡 状 流 一 一 弹 

状 流 转换 准则 , 指出 气泡 直径 大 小 Ds 影响 了 临界 空 
泡 份 额 的 值 : 


j CD j 2 (9882 jo) 0) 
aC, o f 


其 中 ， a 20.2245 < D,» 


S 
a= +0. 1%D, xs 
20D, : E 


a-0.353» 3p 


Taitel (1980) (和 Sadatomi (1982) 上 相继 提 
出 了 适用 范围 较为 宽泛 的 转换 关系 式 ， 形 式 上 分 别 
为 
j 2 0.56, —1.15CogAp / p" (Taite) (3) 


. 1 . 0.35 , 
h= f53c, Ja- C, J 2D, (Sadatomi) (4) 
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可 以 看 出 ， 基 于 漂移 流 模型 提出 的 转换 准则 大 
部 分 含有 漂移 流 分 布 系数 Co 因此 需要 得 到 不 同 矩 
形 通 道内 分 布 系数 的 值 。 
3. 2. 2 矩形 通道 内 分 布 系数 的 计算 
漂移 流 的 形式 一 般 为 : 


UM (0)= 6 U)* Vs (5) 


其 中 <je> 为 平均 气相 折算 速度 ，< a > 为 平均 空 
泡 份额 , < 产 为 折算 速度 ,Vs 为 平均 漂移 速度 。 Zuber 
和 Findley[9] 较 早 提出 了 关于 分 布 系数 的 计算 形式 


C, =1.2-0.2\JD,/D ， 随 后 Jones 和 Zuber', 


IshiiG 分 别 对 于 圆 管 和 窄 通道 提出 了 分 布 系数 的 计 
算 形式 : Co=1.2 and Zuber ) 和 


C, =1.35-0.35,[p,/ p, (shii) 


( Jones 


实验 中 得 到 的 实验 段 的 总 压 降 A P 由 重 位 压 
降 、 摩 擦 压 降 和 加 速 压 降 组 成 。 由 于 通道 是 等 截面 
的 ， 在 绝热 条 件 下 ， 加 速 压 降 可 以 忽略 ， 摩 擦 压 降 
A PP 经 由 Chisholm C1 方法 计算 得 到 ， 则 重 位 压 降 
为 


AP, - AP-AP, - 


而 流 道内 的 截面 含 气 率 已 知 时 , 重 位 压 降 亦 即 : 


AP =| p(1- (a) p, (o) | eL (7) 
其 中 ，< a > 为 两 测 压 点 间 平 均 空 泡 份 额 ,LL 为 
两 测 压 点 间距 ， 于 是 利用 压 降 法 计算 出 的 平均 截面 


A UE 


AP - AF. 
sL (8) 
Pı Pe 


PL 
(a)= 


实验 得 到 的 结果 如 图 4 所 示 ， 随 着 矩形 通道 间 
隙 的 减少 ， 漂 移 流 模型 的 分 布 系数 有 逐渐 增加 的 趋 
势 ， 在 窄 通 道 范围 ， 分 布 系数 可 以 选用 Ishii 进行 计 
算 ， 当 通道 间隙 大 于 5 mm 时 ， 分 布 系数 可 以 近似 
为 1.2。 


C, =1.35-0.35 fp, / p, 


C =1.2 


C, - 12-0240, / p, 一 


= 通道 一 
通 道 二 
a 通 道 三 


(3) m/s 
图 4 ”和 矩形 通道 分 布 系数 比较 


Fig.4 Comparison of distribution parameter in rectangular 


channels 

3. 2. 3 泡 状 流 一 一 弹 状 流转 换 准则 评价 

根据 得 到 分 布 系数 的 计算 方法 ， 对 四 种 转换 准 
则 进行 计算 ， 误 差分 析 结 果 如 表 1 Hz. TEJVMOR 
形 通道 内 泡 状 流 向 弹 状 流转 换 准 则 的 计算 中 ， 对 尺 
寸 较 小 的 通道 Ishii 模型 的 平均 误差 比较 小 ， 因 此 
Ishii 模型 可 以 用 于 窗 矩 形 通道 ， 而 对 于 较 大 尺寸 的 
矩形 通道 ， 四 种 模型 的 误差 都 较 大 ， 适 用 性 不 高 。 
Ishii 模型 通过 比较 气泡 尺寸 与 矩形 管道 间隙 得 到 变 
化 的 空 泡 份额 ， 而 其 余 四 种 模型 则 认为 临界 空 泡 份 
额 保持 不 变 ， 因 此 临界 空 泡 份额 的 计算 是 泡 状 流 向 
弹 状 流转 换 准则 的 关键 。 

表 1 模型 平均 绝对 误差 及 平均 误差 比较 (单位 : €) 


Table 1 Mean deviation and average deviation 


1 


性 


calculated for different correlations (unit: %) 
通道 尺寸 (mm) Ishii Jones-Zuber Taitel Sadatomi 
s=0. 3, w=12 24.7 29.6 50.3 53.4 
s=0. 6, w=12 66.6 70.8 81.5 83.1 
s=], w=12 52.9 65.3 66.2 69.3 
s=1.6, w=43 10.3 13.5 61.2 66 
s=2.4, w=40 21 46.5 29.9 33.6 
s=3, w=43 43.6 18.7 50 56.4 
s=5, w=40 64.1 59.4 79.6 95 
s=10, w=43 62.4 84.9 99 100 


3. 3 泡 状 流 一 一 弹 状 流转 换 准 则 修正 

泡 状 流向 弹 状 流转 换 的 原因 ， 经 典 理论 认为 是 
气泡 间 的 相互 磁 撞 、 合 并 导致 气泡 的 尺寸 增加 ， 以 
至 于 超过 通道 的 直径 而 形成 气 弹 的 结果 ， 因 此 气泡 
的 初始 尺寸 对 于 流 型 转换 的 临界 空 泡 份额 是 有 很 大 
HJE, Ishii 模型 在 奉 抢 形 通道 中 的 适用 性 很 
好 ， 这 主要 原因 是 在 窜 和 矩形 通道 的 实验 研究 中 ， 气 
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泡 初 始 尺 寸 基本 一 致 ,对 于 流 型 转变 容易 形成 共识 ; 
而 在 较 大 间隙 的 通道 ， 由 于 气 液 混合 方式 不 一 致 ， 
造成 气泡 初始 尺寸 相差 较 大 ， 模 型 一 认为 临界 空 泡 
份额 仅 在 0.2-0.3 变化 , 因此 有 必要 在 较 大 间隙 的 和 珑 
形 通道 范围 对 Ishii 模型 的 临界 空 泡 份额 计算 式 进 
行 修正 

文献 [12] 通 过 对 气泡 碰撞 对 于 泡 一 一 弹 状 流转 
换 进 行 数值 计算 ， 得 到 了 一 条 通用 的 、 关 于 气泡 初 
台 尺 寸 影响 气泡 碰撞 率 的 通用 曲线 ， 研 究 表明 当 n 
(7=DwD) 在 0.04-0.4 范围 内 , 临界 空 泡 份额 与 气 
泡 初始 尺寸 的 关系 是 


æ —1.18[(0.664 -- 0.8767) ' — 1] (9) 


将 通过 高 速 摄影 捕捉 到 通道 二 与 通道 三 的 气泡 
初始 尺寸 应 用 到 公式 〈9) 中 ， 并 用 Ishi 的 转换 准 
则 关系 式 的 公式 进行 计算 ， 得 到 新 的 泡 一 一 弹 状 流 
转换 边界 ， 与 实验 值 的 对 比 情况 如 图 5 所 示 ， 实 验 
结果 与 预测 结果 误差 很 小 ， 因 此 在 较 大 间隙 的 矩形 
通道 中 可 以 采用 公式 〈9) 计算 临界 空 泡 份额 。 


结果 
首 三 预 测 结果 


0.1 


气相 折算 速度 (m/s) 
图 5 实验 值 与 预测 值 比较 


Fig. 5 Comparison between experimental data and new 


correction 
4 结论 
CD 基于 漂移 流 思 想 , 对 三 个 实验 段 的 分 布 系 
数 进行 计算 ， 结 果 表 明 在 罕 矩 形 通 道中 可 以 选用 
Ishii 提出 的 分 布 系数 计算 式 ， 而 常规 矩形 通道 中 分 
布 系数 为 1.2。 
(2) 矩形 通道 间隙 大 于 2 mm 时 ， 随 着 通道 水 


力 直径 的 增 大 ， 转 换 边 界 左 移 。 
(3) 在 罕 通 道 范围 可 以 采用 Ishii 提出 的 漂移 
流 公式 作为 流 型 转换 边界 。 
(A) 常规 矩形 通道 中 , 转换 准则 中 的 临界 空 泡 


» 


Tar 


份额 与 气泡 初始 尺寸 有 关 ， 使 用 赵 建 福 计算 临 界 名 
泡 份 额 的 计算 式 修正 Ishii 计算 空 泡 份额 的 公式 , 转 
换 边 界 的 预测 与 实验 结果 符合 程度 较 好 。 
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